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1 .   INTRODUCTION 

Kato"*"  has  demonstrated  the  possibility  of  determining  upper  and 
lower  bounds  on  cot  r\,    where  t\   is  the  phase  shift  for  arbitrary  angular 
momentum  for  the  quantum  mechanical  scattering  of  a  particle  by  a  static 
central  potential.   The  basic  interest  of  the  method  derives  from  the 
fact  that  the  expression  for  either  bound  includes  a  variational  trial 
function,  and  gives  an  error  which  is  of  second  order  in  the  difference 
between  the  trial  function  and  the  exact  function.   The  bounds  are 
therefore  also  variational  estimates,  and  might  be  termed  variational 
bounds.   In  subsequent  studies, ^'^^    the  case  of  angular  momenta 
greater  than  zero  was  considered  more  intensively,  some  slight  exten- 
sions and  improvements  in  the  formalism  were  effected,  and  some  pre- 
liminary results  were  obtained  for  the  applicability  of  the  method  to 
scattering  by  compound  systems  with  emphasis  on  the  low  energy  scatter- 
ing of  positrons  by  Hydrogen  atoms. 

Static  central  potential  scattering  problems  involve  an  infinite 
ntimber  of  channels  which  are  however  uncoupled,  a  consequence  of  con- 
servation of  angular  momentum.   In  other  words,  an  expansion  Into  partial 
waves  leads  to  a  dlagonalization  of  the  scattering  matrix  S 
before  the  detailed  analysis  is  begun.   Furthermore,  each  an^lar  momen- 
tum channel  wave  function  is  completely  characterised  by  one  real  number, 
the  real  phase  shift  for  that  angular  momentum.   More  generally,  the 
Kato  method  can  be  applied  whenever  the  S  matrix  can  be  diagonallzed  on 
some  a  priori  basis  (generally  this  will  be  through  symmetry  considera- 
tions); the  resultant  elgenmodes  can  then  each  be  characterized  by  one   • 


"5 
real  number.   Examples   include  the  one  dimensional  scattering  problem 

for  an  even  potential^  V(x)  =  V(-x),  and  the  problem  of  electromagnetic 
waves  in  a  rectangular  wave  guide  which  has  an  obstacle  symmetric  with 
respect  to  some  plane  perpendicular  to  the  axis  of  symmetry.   The  latter 
case  involves  vector  wave  functions  of  three  independent  variables  rather 
than  scalar  wave  functions  of  one  independent  variable,  but  this  causes 
no  essential  changes.   In  both  cases,  the  eigenmodes  are  the  even  and 
odd  standing  wave  solutions. 

In  this  paper,  we  will  consider  multi- channel  scattering  in  which 
the  scattering  matrix  cannot  be  dlagonalized  a  priori .   More  foimally, 
we  are  concerned  with  elastic  multi-channel  scattering  processes.   We 
shall  cast  the  Schroedinger  equation  into  a  form  in  which  the  potential 
is  a  real  symmetric  N  X  N  matrix  and  the  wave  function  is  an  N  X  1  matrix, 
though  it  is  not  necessary  to  do  so.   The  equation  for  neutron  proton 
scattering  with  tensor  forces  is  conveniently  written  in  matrix  form 
with  N  =  2,  and  we  shall  see  in  Section  2  that  the  same  is  true  for  the 
problem  of  one -dimensional  scattering.   It  is  our  purpose  to  determine 
upper  and  lower  bounds  on  the  elements  of  the  S  matrix,  the  N  X  N  ma- 
trix which  relates  the  N  outgoing  waves  to  the  N  incoming  waves.   The 

2 
elements  of  S  are  in  general  complex,  so  that  there  are  2N  real  quan- 
tities in  all,  but  it  is  well  known  that  only  ^N(N+1)  of  these  real 

7 
quantities  are  independent. 

The  procedure,  in  analogy  with  that  of  Kato,  is  to  derive  an  inte- 
gral variational  principle  which  is  in  fact  an  identity  and  then  to 
bound  the  explicitly  exhibited  second  order  term.  Upper  and  lower  bounds 


are  first  obtained  on  the  elements  B. .  of  B,  where  the  susceptance  matrix 
B  is  the  inverse  of  the  Cayley  transform  of  S.   One  can  then  also  obtain 
bounds  on  the  elements  of  S,  or  on  still  another  set  of  ^N(N+1)  quantities 
which  characterize  the  scattering,  the  N  eigenphase  shifts  and  the  s^CNyl) 
mixing  parameters . 

Numerical  examples  are  calculated  in  Section  6  for  a  one -dimensional 
attractive  square  well  asymmetric  with  respect  to  the  origin. 

Because  of  the  enoiinous  advantages  of  variational  bounds  over  varia- 
tional estimates,  it  is  of  great  interest  to  study  the  limits  to  which 
the  Kato  method  can  be  extended.   There  are  a  number  of  obvious  directions 
in  which  to  proceed.  We  will  list  some  systems  of  interest,  starting  with 
those  already  considered. 

(1)  Systems  which  involve  the  scattering  of  a  particle  by  a  static 
central  potential.  We  are  then  of  course  concerned  with  one  channel  pro- 
cesses. 

(2)  Systems  which  for  sufficiently  low  energy  of  the  incident  particle 
involve,  as  does  (l),  only  one  (elastic)  channel,  or  a  priori  sepajrable 
(elastic)  channels,  but  which  are  compound.  We  will  refer  to  a  system  as 
compound  if  it  has  spatial  (as  opposed  to  finite  dimensional  spin)  inter- 
nal degrees  of  freedom. 

(5)   Systems  which  involve  the  scattering  of  a  particle  by  a  potential 
which  is  static  (the  scattering  system  has  no  spatial  degrees  of  freedom) 
and  which  can  be  represented  by  a  hermitian  N  X  N  matrix.   It  is  assiimed 
that  the  N  >  1  elastic  channels  cannot  be  imcoupled  in  any  a  priori  way. 


{k)      Systems  which  are  compound  and  vhich  involve  N  >  1  elastic 
channels . 

(5)   Systems  which  involve  inelastic  scattering. 

From  a  purely  formal  point  of  view^  case  (2)  is  included  in  the  ori- 
ginal paper  of  Kato.    Unfortunately,  the  application  of  the  method  re- 
quires the  evaluation  of  lower  bounds  on  two  eigenvalues,  a  and  p,  which 
arise  in  an  associated  eigenvalue  problem.   While  these  eigenvalue  bounds 
are  relatively  simple  to  obtain  for  scattering  by  a  static  central  poten- 
tial, no  methods  for  obtaining  bounds  on  a  for  scattering  by  a  compound 
system  have  been  developed. 

The  determination  of  a  lower  bound  on  p  gives  an  upper  bound  on 
(-k  cot  T\)       .      The  determination  of  this  upper  bound  through  the  Kato  ap- 
proach has  however  now  been  rendered  largely  academic  by  the  development 
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of  a  technique  '      which  does  not  require  the  evaluation  of  matrix  elements 

of  the  square  of  the  Hamlltonian,  as  does  the  Kato  approach.   This  tech- 
nique, incidentally,  is  a  fortiori  applicable  to  (l). 

The  Kato  method  as  originally  formulated  was  explicitly  restricted 
to  one  channel  processes,  and  could  therefore  not  be  applied  to  (3)^  (^)^ 
or  (5).   The  present  paper  concerns  itself  with  the  development  of  exten- 
sions of  the  Kato  method  which  will  allow  it  to  be  applied  to  (5);  further, 
explicit  methods  will  be  given  which  will  often  allow  one  to  obtain  lower 
bounds  on  a  and  on  (3.   [These  extensions  can  also  be  used,  in  principle, 
for  (k-) ,    and  can  probably  be  extended  to  (5)^  but  these  cases  have  not  as 
yet  been  seriously  examined."]   Now  work  in  progress  shows   that  the  tech- 
niques referred  to  can  be  applied  to  {3),    {^) ,    and  (5).   One  therefore  has 


a.  very  much  simpler  method  for  obtaining  an  upper  bound  on  (-k  cot  T])" 
for  case  (3)  than  can  be  realized  by  the  present  paper.   The  present  work 
is  still  justified  because  one  would  of  course  like  very  much  to  be  able 
to  detennine  the  other  bound.   (it  is  of  course  also  preferable  not  to 
have  to  truncate  the  potential,)   Since  it  entails  very  little  additional 
labor,  we  will  actually  obtain  both  bounds  for  (3)- 

2.   MATRIX  FORMULATION  OF  THE  VAEIATIORAL  PRINCIPLE 
We  begin  by  showing  how  the  one  dimensional  scattering  problem  with 
potential  V(x)  can  be  written  as  a  differential  equation  in  matrix  form. 
We  have 


2  2 
-ft  d  u(x) 

2m   ^  2  "^ 
dx 


V(x)  -  E 


u(x)  =  0,  (1) 


where  we  take  V(x)  to  be  real  and  to  vanish  more  rapidly  than  l/x  for 
|x|  -►  00.   It  is  evident  that  there  are  two  channels,  corresponding  to 
X  >  0  and  to  x  <  0,  respectively,  which  are  coupled  by  the  potential, 
and  consequently  tliree  independent  parameters  are  required  to  specify  the 
asymptotic  effects  of  the  scattering  process. 

It  is  convenient  to  recast  Eq.  (l)  in  matrix  form  so  that  the  three 
unknown  parameters  appear  as  the  three  independent  elements  of  a  single 
matrix  quantity.   This  can  be  accomplished  by  expressing  u(x)  and  V(x) 
as  sums  of  even  and  odd  functions. 


u(x)  =  u^(x)  -  u^(x),    V(x)  =  Vg(x)  +  V^(x),  (2) 


vhere 


u   (x) 


i[u(x)  +u(-x)J 


(3) 


and 


V   (x) 


=  iLv(^)  ±  V(-x)] 


(M 


V  and  V  will  clearly  also  vanish  more  rapidly  than  l/x  as  |x| 


If 


we  replace  x  "by  -x  in  Eq.  (l)^  add  and  suhtract  the  derived  equation  and 
Eq.  (l)^  and  use  Eqs .  (5)  and  {^) ,    we  obtain  the  single  matrix  equation 


(H-E)u  =  - 


2m 


2  d^u 


dx 


u  =  0 


(5) 


where 


u  = 


u 


V  = 


V 


:6) 


Note  that  the  channels^  which  are  still  coupled  in  general,  now  correspond 
to  the  even  and  odd  portions  of  the  wave  function  rather  than  to  x  >  0  and 
X  <  0.   We  can  restrict  the  range  of  the  independent  variable  to  x  ^  0  since 
the  components  of  u  are  even  and  odd  functions  of  the  coordinate.   The  bo\m- 
dary  conditions  at  the  origin  are,  obviously, 


u  =  0  ,   du  /dx  =  0 
o      '  e' 


(7) 


We  see  then  that  one  can  conveniently  formulate  one  dimensional 
scattering  in  terms  of  a  Schroedinger  equation  vlth  V(x)  a  real  her- 
mitian  2x2  matrix.   For  the  formal  development  vhlch  follows,  we  will 
not  restrict  ourselves  to  matrix  potentials  of  the  above  form  but  rather 
we  will  take  V(x)  to  be  an  arbitrary  real  hermitian  2x2  matrix  whose 
elements  fall  off  faster  than  l/x.   The  indices  will  therefore  be  1  and 
2  rather  than  e  and  0.   The  results  to  be  obtained  are  actually  valid 
for  V  an  N  X  N  matrix,  that  is,  for  an  N  channel  process;  we  set  N  =  2 
only  as  a  matter  of  convenience  since  it  slightly  simplifies  the  discus- 
sion and  the  notation.   Many  elastic  scattering  processes  involving  com- 
pound systems  can  be  written  in  a  slightly  generalized  form  of  Eq.  (5)- 

We  will  now  derive  a  variational  principle  for  scattering  which  can 
be  described  by  Eq.  (5)-   The  asymptotic  form  of  u(x)  can  be  written  as 

.  ,      ^-i(kx+G)       +i(kx+e)  (8a) 

k^  .   2mE/*2.  (8b) 

The  utility  of  the  parameter  G,  where  0  S  0  <  «  but  is  otherwise  arbitrary, 
will  become  apparent  later.   It  is  convenient  to  rewrite  Eq.  (8a)  in  a  form 
that  involves  standing  waves  rather  than  traveling  waves,  since  we  will 
then  be  able  to  deal  with  hermitian  rather  than  unitary  matrices. 


ufx)  ->  e^  cos(kx  +  e)  -  ij„  sin(kx  +  ©),  (9a) 


The  amplitude  vector  for  the  outgoing  waves,  h   Is  related  to  that  for 
the  Incoming  vaves_,  a  _,  by  the  scattering  matrix  S^  through 

The  vectors  e   and  j   are  related  hy  the  susceptance  matrix  B   (so 
called  because  it  is  analogous  to  the  susceptance  matrix  in  network  theory), 


From  Eq.  (9)  and  (lO)  it  can  be  seen  that  B   is  the  inverse  of  the  Cayley 

**«y 

transform  of  S  ^ 

Bq  =  i(l  +  Sq)(1  -  Sq)-^  .  (12) 

12 

It  can  be  shown  that  S   Is  unitary  and  symmetric   and  as  a  consequence 

*A^  y 

that  B  is  hermltia,n  and  symmetric  and  therefore  real.   Equation  (ll)  can 
-"•9 

now  be  used  to  eliminate  j   in  Eq.  (9a),  giving  rise  to 

u(x)  -*  e^  cos(kx  +  9)  -  B^  e^  sin(kx  +  9).  (13) 

«.^     ~»y  ^-9  -^y 

We  now  Introduce  a  trial  wave  function  u  (x)  which  is  required  to 
have  the  asymptotic  form 

u_^(x)  ->  e^  cos  (kx  +  9)  -  B^,  e   sln(kx  +  9)  (l^) 

M\X,     ^~*y  -^yc  *** y 

as  X  ^  +00,  that  is,  the  amplitude  vector^  is  the  same  as  that  of  the 
exact  function,  but  the  unknown  B^  is  replaced  by  B^^ .  In  addition,  we 
impose  the  boundary  condition 

u  du  /ix  =  (du  /dx)u 
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at  X  =  0,  and  we  require  that  u  and  du  /dx  be  continuous  in  the  range 

«*t      '^t 

0  <  X  <  00  .   It  is  to  be  noted  that  all  of  the  boundary  conditions  im- 
posed on  u  are  also  satisfied  by  the  exact  function.   In  order  to  arrive 
at  a  variational  principle^  consider  the  relation  (all  integrals  extend 
from  0  to  co^  and  the  limits  will  henceforth  be  omitted) 


■(2m/-«^) 


utn  u  -  (H  u)'l"u 


dx  = 


.f± 


*^  dx  ;>-t   dx  ~  -t 


(15) 


Using  Eq.s .    (13)    and   (5)^    the  boundary   condition   at   the   origin^    and  the 
symmetry  of  B^    we   obtain   from  Eq.    (l5) 


£eV~e  =  s/?0t^0-^^/*^)    /    ^t^^l-'^^^t'^ 


+   (2m/«^0       /     v"''(H-E)w  dx, 


w   =  u      -   u 

*^  ^t  — 


(l6a) 
(l6b) 


The  first  two  terms  on  the  right  hand  side  of  Eq.  (l6a)  can  be  evalu- 
ated directly  from  the  trial  function,  and  the  third  term  is  of  second 
order  in  the  difference  between  the  trial  function  and  the  exact  function. 
The  first  two  terms  therefore  constitute  a  variational  approximation  for 
e  '''b  e  .   However,  the  approximation  Is  in  general  neither  an  upper  nor  a 
lower  bound  on  the  exact  value.   In  the  next  section  we  will  derive  rigor- 
ous second-order  upper  and  lower  bounds  on  the  error  term,  the  third  term 
on  the  right  hand  side  of  Eq.  (l6a). 


11 


5.   RIGOROUS  BOUNDS  ON  THE  ERROR  TERM 
In  this  section  we  will  derive  rigorous  upper  and  lower  bounds  on 
the  error  term  in  Eq.  (l6a),  following  a  procedure  analogous  to  that  de- 
rived by  Kato.   We  begin  by  considering  the  associated  eigenvalue  prob- 
lem. 


(H-E)sz$(x)-  (ip(x)jZ^(x)  =  0. 


(IT) 


In  Eq.  (17),  P(x)  is  required  to  be  a  real,  positive  definite  hermltlan 
matrix  each  of  whose  elements  must  fall  off  faster  than  l/x  as  |x|  -^-  0°  . 
In  particular,  we  are  interested  in  discrete  eigenvalues  |a  ,  of  which 
there  will  be  an  infinite  number,  defined  by  the  condition  that  H  is  a 
hermltian  operator  when  operating  on  the  corresponding  eigenf unctions, 
^  ,  or  on  w   [defined  by  Eq.  (l6b)2,  that  is. 


^^  E  ^   -  {E  ^  y  ^ 

Zm     «  J:n    ~  ~m      „n 


dx  =  0 


(l8a) 


w"*"  H  (Z^   -  (H  w)"!"  (Z^ 
%«  *«  n    -^   ***   J,  n 


dx  =  0  . 


(18b) 


It  follows  from  the  hermltlan  character  of  H  that  the  associated  eigen- 
functions  are  orthogonal  with  respect  to  the  weight  factor  p  and  that 
the  associated  eigenvalues  are  real. 


^J   p  (Z^^  dx  =  {^\/2m)b. 


.m  ^  -.n 


mn 


(19s 


•^n    ^n 


(19b) 
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The   normalization  in  Eq.    (l9a)    has  been    chosen  for   later   convenience. 
In   addition  we   require   that   the   ^      satisfy  the  boundary  conditions   at 


X  =   0, 


i^  ^L   I   ^  =    ^^Q  I  ^)   in   '  ^20a) 

(Z^tdw/iiczz    (djz^t/  dx)w    .  (20b) 


_.n 


In  analogy  vith  Eq.  (13)^  '^     has  the  asymptotic  form  as  x  ->■  oo  , 


.n 


<z^  -»  e   cos  kx  -  Bfu  )  e   sin  kx  .  (21) 

^n   «,n         ~^^n  -,n  ^ 

From  Eqs .  (ij)^  (l^)  and  (l6b)  it  follows  that  the  asymptotic  form  of  w 
is 

w-  (Bq  -  Bq^)  e^  sin(kx  +  9).  (22) 

Using  Green's   theorem  and  Eqs.    (20),    (21),  and   (22)    to  evaluate   the 
integrals   in  Eqs.    (l8)    we   find  that   Eq.    (l8)    can  be   satisfied  only  if 
B(|j.   )    satisfies 

B(Li    )    e     =   -    cot  e  e      .  (25) 

M,^'^n     -n  v.n  ^       ' 

Thus  the  associated  eigenfunctions  are  also  eigenmodes,  i.e.,  their  ampli- 
tude vectors  e   are  eigenvectors  of  the  susceptance  matrix  B(|a  ).   The 
A*n  "-^  n 

\i     are  just  the  values  of  [x   for  which  one  or  another  of  the  eigenvalues 
of  the  susceptance  matrix  is  equal  to  -  cot  9. 

We  will  now  proceed  to  derive  a  monotonicity  theorem  required  in  sub- 
sequent developments.   In  principle,  we  can  find  solutions  i^(x)  of 


15 


Eq.  (17)  for  arbitrary  (j.  if  we  no  longer  impose  the  boundary  conditions 
(18).  In  particular^  we  can  find  2  eigenmodes  ^.(x)  (i  =  1,2)  for  an 
arbitrary  value  of  |j  which  are  distinguished  by  their  asymptotic  forms, 


(n)  -  eAu) 


<vi 


cos  kx  -  ^.(ji)  sin  kx 


(24; 


,th 


In  Eq.  (2^)  ^.(n)  is  the  1   eigenvalue  of  B(|_i).   It  should  be  noted 
that  the  eigenmode  solutions  of  Eq.  (l?)  are  employed  here  only  in  the 
formal  development.   Presumably,  it  is  more  difficult  to  find  such  solu- 
tions than  to  solve  the  original  problem,  Eq.  (5)-   At  this  point  it  is 
convenient  to  relabel  the  associated  eigenfunctions  with  two  indices. 

The  associated  eigenvalues  u.  .  are  then  defined  by  the  relation 

ni 


f. (n  . )  =  -  cot(9  +  nn) 
^1  ni 


(25) 


where  the  index  i  labels  the  eigenmodes  and  the  corresponding  associated 
eigenfunctions  have  the  asymptotic  form 


^ni   «,i  nx 


cos  kx  +  cot(0  +  nn)  sin  kx 


(26; 


Consider  the  eigenmode  solutions  of  Eq.  (17)  for  two  different  values 
of  |i,  to  be  denoted  by  ja'  and  |a",  without  applying  the  boundary  conditions 
(18),  and  form  the  expression 


•(2m/-fi^) 


.ut, 


(^  +(^i•)H  ^.iu-")    -   (H  /.(^'))7.(^^" 


dx 


{2m/^%'*   -  (i")   /  ^i■^(^')  P  iji^")^^ 


li 


^"^   dx  ~j       dx  «i  ^    „j 


(27) 
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■where  we  have  made  use  of  Eqs .  (5)  and  (l?).   Using  Eqs .  (20)  and  the 
asymptotic  form  of  (i^.(i-i)  given  by  Eq..  (24),  we  may  rewrite  Eq.  (27)  as 


(2m/*^)(n'*  -  n")  f    s2$.+(n')  P  ^Ali")^ 


=  k  e  +(n>)  e  (n") 

•iMl  -»J 


q(n')    -    ^j(^i") 


(28) 


Consider  Eq.  (28)  with  ji'  and  |j."  real  but  otherwise  arbitrary  and 
with  j  =  1.   Let  the  eigenphase  shift  5.(|j.)  be  defined  by 

^.(n)  =  -  cot  8.(n).  (29) 

Substituting  Eq.    (29)    in  Eq.    (28)    and  taking  the   limit   |j."   -*  n',    we 
obtain 


db./dii  =   (2m/«^)sln^S.       /      ^"''   P   jZ^.dx  /    (ke.+e.)    g  0.  (50) 


.2,     .    2, 


The  Inequality  in  Eq.  (50)  follows  from  the  positive  definite  character 

of  p  .   Thus  the  eigenphase  shifts  are  monotonlc  increasing  functions 

of  [1  (for  k  ^  0). 

We  now  proceed  to   derive   rigorous  bounds    on  the   error  term  in  Eq. 

(16a).      Suppose  that  f  and  g  are   functions  which  either  have  the   same 

phase   shift   as   the   ^    .    or  fall  off  asymptotically  faster  than  l/x.      We 

define  expansion  coefficients   a   .    and  b    .    by  the  relations 
■^  mm 


a    .    =    (2m/«^)      /      ^    .■''   P   f  dj-   ,  (51a) 

m        ^      '         '  Jim     ^  '-         ' 
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b  .  =  (2in/ft^)   /   (Z^  ."*"  p  g  dx  .  (Jib) 


ni    '         /  J.  ni 


It  can  then  be  shown  that  the  associated  eigenf unctions  form  a  complete 

13 

set    in  the  sense  that 


f "*"  p  f  dx  =  (ft\/2ra)    ^   a  .*  a  .  (32a) 


»  ,„  -  :-    ni   ni 

n,i 


f  %  g  dx  =  (^V2m)   ^   a  *  b  (32b) 

n,i 


g"^  p  g  dx  =  (ftVsm)   ^   b   *  b  (52c) 


ni   ni 
n,l 


If  we  now  set  f  =  w  and  g  =  -p   (H-E)w  ,    it  follows  from  Eqs .  (31)  that 


^       >U4       /K 


a  .  and  b  .  are  related  by 
ni      ni 


b  .  =  -(2m/fi^)   /  <i>    .''"(H-E)w  dx  =  -(2mM\)   /  Rh-E))!^  .1  "^w  dx 


-\i   .(2m/*^)   /   (Z^  ."''  p  w  dx  =  -u  .  a  .,  (33) 

^ni  ^   '      /  JT ni  -«.  *i-       ni  ni 


•k  j-i-^  ■«« 


provided  p  is  chosen  such  that  g  falls  off  asymptotically  faster  than  l/x. 
Eqs .  (52)  can  now  be  written  in  the  form 
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w+  (H-E)w  dx  =  (*^/2m)   ^   u  "^  b   *  b    ,  (5^b) 

A«     AM     AW  n  n 


r  [(H-E)w]'^  p"^(H-E)w  (be  =  (*V2m)   ^   V^^  V  *       ^^""^ 


Let  a  be  the  smallest  positive  eigenvalue  and  -P  be  the  smallest 
9  " 

(in  absolute  value)  negative  eigenvalue.   It  follows  from  Eqs .  (5^)  that 


-p  "^  e/  S   r  w"^(H-E)w  dx  S  a  "^  c/  ,  (35a) 


e^=  [(H-E)w]^  p  ^(H-E)w  dx  .  (55b) 

9         ,  /         ^^      /»«.      v»      *.»     A< 


We 


have  thus  obtained  rigorous  upper  and  lower  bounds  on  the  error 


2 
term  in  Eq.  (l6a).   The  quantity  €q  is  of  second  order  in  the  difference 

between  the  trial  function  and  the  exact  function,  and  can  be  calculated 

from  the  trial  function  since  it  follows  from  Eqs .  (5)  and  (l6b)  that 


(H-E)w  =  (H-E)u   .  (36) 
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The  problem  still  remains  of  course  of  obtaining  lower  bounds  on  a  and 
on  p  .  It  is  to  be  recalled  that  it  is  possible  to  obtain  close  bounds 
on    /  "w  (H-E)wdjt   even  with  quite  crude  lower  bounds  on  a^  and  p  , 

J  U^  AM  fM^  9  0 

2 
since  e    is  a  second  order  quantity. 


h.      BOUNDS  ON  a  AND  ON  p^ 
In  this  section  we  will  indicate  a  rather  general  method  for  ob- 
taining lower  bounds  on  a  and  on  p   for  cases  covered  by  Eq.  (5)-   The 

y       y 

monotonicity  theorem,  Eq.  (50),  states  that  the  eigenphase  shifts  5.((j.) 
are  monotonic  increasing  functions  of  \i.      The  labeling  of  the  eigenphase 
shifts  is  essentially  arbitrary,  but  the  application  of  this  theorem 
requires  that  we  keep  track  of  the  labeling  as  \x   is  varied.   Since  [x   is 
3.   continuous  variable,  this  presents  no  difficulty  except  at  crossing 
points  where  5.(|a)  =  6.(ti),  i  ^   j.   However,  we  are  interested  in  the 
weaker  statement 


< 


^n-1,1    ni    n+1,1 


(57) 


where  \i    .    is  the  value  of  \x   for  which  5.(|a  .  )=  0  +  njt .   Eq.  (57)  holds  in 

spite  of  the  ambiguity  in  labeling  at  crossing  points,  since  no  matter 

which  curve  is  continued  into  which,  both  continue  to  rise  raonotonically. 

The  relationship  between  a  and  p   and  the  eigenphase  shifts  is  illus- 

y      y 

trated  for  a  typical  case  in  Figure  1.  If  we  define  a  and  -p  to  be 
the  smallest  positive  eigenvalue  and  smallest  (in  absolute  value)  nega- 
tive eigenvalue,  respectively,  associated  with  the  i    eigenmode,  then 

a  is  the  smaller  of  the  a  .  and  p^  the  smaller  of  the  p   .   Thus  there 
0  01      9  01 

is  only  one  a  and  one  p  . 
9  0 
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S.(/i) 


8.(/i)    Sgl/i) 


d-TT 


Figure  1 
Relationship  between  a  and  p   and  the  eigenphase  shifts 
of  the  associated  eigenvalue  problem,  5.(|a).   The  eigen- 
phase shifts  of  the  original  problem  are  t\     and  j]^. 
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Figure  1  demonstrates  the  relationship  between  a  and  (3   and  the 
eigenphase  shifts,  but  these  eigenphase  shifts  are  presumably  not  known. 
We  must  therefore  devise  a  method  for  obtaining  lower  bounds  on  a  and 
P  .   To  this  end,  we  introduce  a  comparison  potential  V  such  that 

y  *«  c 

A  V  =  V   -  V  is  non-positive  and  for  which  the  associated  eigenvalue 
problem  can  be  solved  exactly,  or  for  which  bounds  can  be  obtained.   We 
will  now  show  that  the  eigenphase  shifts  corresponding  to  the  comparison 
potential  are  upper  boimds  on  the  eigenphase  shifts  corresponding  to  the 
true  potential  for  any  given  value  of  |a.  Under  the  substitutions 

H  — >  H  -  la  P  (  38a) 

*^X         Av        A^ 

H  -  n  P  — >  H-nP+XAV,  (38b) 

M„  AAA  AM  'Vb'  VU^     -Ua 

the  derivation  which  led  to  Eq.  (30)  gives 

d6./dX  =   -(2m/R^)sin^5.       /      ^  "^  A  V  ^    dx  /   (ke.  "*"  e.  )   g  0    .  (39) 

1  X         !         ftui        Om    au    mi.  ***  i        •*A  1 

The  inequality  follows  from  the  non-positive  character  of  A  V.   The  eigen- 

AAA.  *A, 

phase  shift  5  .{\i)    corresponding  to  the  comparison  potential  is  obtained 
when  X  =  1,  and  5.(^)  when  X  =  0.   It  follows  from  Eq.  (39)  that  5  .(^) 
is  an  upper  bound  on  6.(n), 

8^i(n)  ^  6.(^).  (ho). 

As  before,  we  are  required  to  keep  track  of  the  eigenphases  as  X  is  varied 
continuously  from  zero  to  one.   The  only  difficulty  occurs  at  a  possible 
crossing  point,  but  B.(|a,X)  =  5.(|a,X)  already  implies  that  these  eigen- 
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phase  shifts  axe  greater  than  the  true  values  and  they  can  only  increase 
as  X  approaches  one.   Eq.  (^O)  is  therefore  valid  in  any  case. 

By  introducing  another  comparison  potential  such  that  A  V  is  non- 
negative^  it  can  be  shown  in  a  similar  fashion  that  the  corresponding 
eigenphase  shifts  6   (|a)  are  lower  boimds  on  the  5.((i)^ 

B^.(u)  §  5.(^)  .  (J+1) 

The  way  in  which  the  bounding  curves  obtained  from  the  comparison 

potentials  can  be  used  to  calculate  lower  bounds  a   and  p   on  a  and 

on  p  ,  respectively,  is  illustrated  for  a  typical  case  in  Figure  2.   It 
y 

is  clear  from  Figure  2  that  if  the  comparison  potentials  do  not  approxi- 
mate the  actual  potential  fairly  closely,  it  may  not  be  possible  to  ob- 
tain lower  bounds  on  a  and  on  p  .   For  example,  if  the  curve  labeled 
6  ,(|j.)  in  Figure  2  were  to  cross  the  dashed  line  5.(n)  =  9  at  a  point 
where  |j.  is  negative,  then  it  would  n  ',  be  possible  to  obtain  a  lower 
bound  on  a  . 

In  general  it  is  possible,  ay      '-   shall  show,  to  calculate  lower 
bounds  on  a  and  on  p  by  using  comj  r:.son  potentials  and  a  weight  func- 
tion  p  which  are  diagonal.   In  that  c  .se,  the  bounding  curves  5  .{\i)    and 
5   (|a)  are  obtained  by  solving  h   singl 2- channel  problems  rather  than  2 

■61 

2-channel  problems.   The  consequences  of  this  result  are  far  reaching, 
in  that  very  few  truly  multichannel  problems  can  be  solved  at  all,  even 
with  simple  comparison  potentials.   On  the  other  hand,  the  possibility 
of  using  diagonal  comparison  potentials  implies  that  it  is  no  more  diffi- 
cult to  obtain  bounds  in  the  multichannel  problem  than  in  the  single- 
channel  problem,  though  it  may  be  much  more  tedious,  and  thus  the  power 
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Si  (/I) 


\S,,WSu2(M' 


Figure   2 


Use   of  comparison  potentials   to   obtain  lower  bounds,    a. 

'     ee 

and  Pq^j,  on  a^  and  (B^,  respectively.   5  .(n)  and  5^.(ia)  are 
derived  from  the  comparison  potentials  and  are  upper  and 
lower  boiinds,  respectively,  on  5.(|a). 
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of  the  method  is  greatly  enhanced.   Furthermore,  the  Kato  method  for 
obtaining  bounds  in  the  single- channel  case  can  be  used  to  obtain  bounds 
on  the  eigenphase  shifts  of  the  comparison  potentials,  or  the  single- 
channel  comparison  potential  problems  could  even  be  solved  numerically. 
We  will  now  consider  the  requirements  that  must  be  imposed  on  the 
elements  of  V  (x)  if  it  is  talcen  to  be  diagonal.   Writing 

('- 

we  have 


^1  -  ^11  -\2 

A  V  =   I  1  (^5) 


^<^    iM 


^21        V^2-  ^22 


It  is  then  trivial  to  show  that  if  A  V  is  to  be  non-negative,  we  must 
have 

^cl  =  ^11  ■"    \\2\         '        \2  ^  ^22  ^    l\2l    '  (^^^ 

while   if  A  V  is   to  be  non-positive,    we  must  have 


^cl  ^  \l   -    l\2l         '        \2  =  ^22   -    1^12!     •  (^5) 


The  above  conditions  are  necessary  but  not  sufficient  for  obtaining  lower 

bounds  on  a  and  on  p  ,  since  we  have  already  seen  that  it  is  possible 
0       y 
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that  no  diagonal  comparison  potential  may  be  sufficiently  close  to  V  to 
be  able  to  determine  bounds. 

5.   METHOD  OF  OBTAINING  BOUNDS 
We  are  now  in  a  position  to  apply  the  foregoing  results  to  the  prob- 
lem of  obtaining  rigorous  upper  and  lower  bounds  on  the  independent  ele- 
ments of  the  scattering  matrix.   From  Eqs.(l6)  and  (55),  we  can  obtain 
expressions  which  are  upper  and  lower  bounds  on  a  quadratic  form  of  the 
susceptance  matrix. 


g  (2me//fi^a.,),  (k6) 


-0  '     91 


[(H-E)u  ]^  p"^(H-E)u,dx  ,  (47) 


where  a   and  (3   are  lower  bounds  on  a  and  on  (3  ,  respectively.   We 
thus  have  variational  bounds  on  e  Be.   The  boimd  property  enables 
us  to  improve  the  approximations  by  incorporating  variational  parameters 
into  u  .  As  opposed  to  the  situation  in  the  usual  variational  calcula- 
tion,  we  do  not  have  to  find  the  stationary  values  of  the  parameters 
(it  is  preferable  but  not  essential)  but  can  rather  evaluate  the  bound 
for  a  few  values  of  the  parameter  and  unambiguously  choose  the  best  re- 
sult.  This  makes  feasible  the  use  of  variational  parameters  which  are 
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not  linear.   0  can  be  thought  of  a  such  a  parameter,  though  in  the  pre- 
sent paper  we  will  not  take  full  advantage  of  this  property. 

As  an  alternative  procedure,  we  could  adjust  the  variational  para- 
meters  to  minimize  the  quantity  e   given  by  Eq.  (^7)-  While  this  proce- 
dure  is  simpler,  in  general  it  yields  inferior  bounds  and  consequently 
it  will  not  be  followed  in  this  paper.   With  reference  to  Figure  2,  9 
should  be  chosen  to  provide  the  greatest  lower  bounds  on  a  and  on  (3  . 
In  the  single- channel  problem  9  can  be  chosen  differently  for  a  and  for 
P  ,  since  the  upper  and  lower  bounds  on  the  phase  shift  are  obtained  by 
solving  two  distinct  problems.   However,  in  the  multi-channel  case  the 

susceptance  matrices  B  corresponding  to  two  different  values  of  9  are 

«^9 

related  in  such  a  complex  manner  that  it  does  not  appear  feasible  to 

obtain  upper  and  lower  bounds  on  the  5  independent  parameters  of  the 

scattering  matrix  by  using  different  values  of  9  for  the  upper  and  lower 

boimds  on  the  quadratic  forms   e  Be. 

*^9  »** 9*'9 

In  order  to  obtain  bounds  on  the  3  distinct  elements  of  B  ,  we  must 

evaluate  the  right  hand  sides  of  Eqs .  (J+6)  with  3  different  choices  of 

the  trial  function  which  differ  in  the  amplitude  vector  e  .   Note  that 

»»9 

the  different  values  of  e  are  not  required  to  be  linearly  independent 

•w9     

(they  of  course  could  not  be)  since  e  appears  quadratically  in  e  B^e  . 

Rather,  the  resulting  inequalities  in  the  3  independent  elements  of  B 

must  be  linearly  independent.   Since  there  are  two  bounds  for  each  choice 

of  e^,  there  are  6  problems  in  all. 
«9 

Since  the  scattering  matrix  S  (the  scattering  matrix  S  as  normally 
defined  is  S^  with  9=0)  can  be  calculated  from  B  ,  the  3  independent 
elements  of  B_  are  siifficient  to  determine  S-   S  can  also  be  calculated 

•"0  t^  AM 
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from  an  alternative  set  of  5  parameters,  the  2  eigenphase  shifts  and  the 
mixing  parameter.   This  alternative  mode  of  description  is  useful  in  many 
applications. 

We  now  examine  the  procedure  for  obtaining  bounds  on  the  eigenphase 
shifts  and  mixing  parameters.   Since  B  is  hermltlan  and  real,  it  can  be 
dlagonalized  by  an  orthogonal  matrix,  T.   The  same  matrix  T  also  dlagon- 
allzes  S,  since  B  is  related  to  S  by  a  rational  fraction.   Because  S  is 
unitary,  its  diagonal  form  can  be  written  as  exp(2iri)  where  t]  is  a  dia- 
gonal  matrix  whose  elements  are  the  eigenphase  shifts.   Note  that  when 
the  eigenphase  shifts  are  defined  in  this  way,  they  do  not  necessarily 
vanish  in  the  absence  of  a  potential  but  rather  assume  values  consistent 
with  the  boundary  conditions  at  the  origin.   The  mixing  parameter  is 
just  the  angle  of  rotation  on  which  T  depends.   It  follows  from  Eqs .  (8a), 
(10)  and  (12)  that  B„  is  related  to  S  by 

B  =i(l^Se-2ie)(i.se-2ie)-l  .  (^8) 

*M  y     <*•   «M      t/ut.      A** 

From  the  foregoing  discussion,  S  can  be  represented  by 

S  =   T  exp(2iT|)f  .  (^9) 

A«*  Ma  *i*        MA 

Therefore,  from  Eq.  (^8)  and  (^9)  B   can  be  written  as 

Mt  y 

B_  =  -T  cot(Ti-0)T  .  (50) 

•"  y   •"    M<   >«• 

The  orthogonal  matrix  T  can  be  written  as 


(51) 


2k 


where  ^   is  the  mixing  parameter.   Solving  Eq.  (50)  and  (51)  for  -i/ ,    t)  , 


and  T\   ,    we  obtain 


^  =  2  cot 


^^922  -  ^eil^  /  (2B^12) 


(52a) 


+  cot 


*(^e22-  ^eil)^°^2t  +  B^^2^in2x^  -^^^^^^   B^^^) 


(52b; 


■q„  =  0  +  cot 


2(Bq22-  BQ^^)cos2t 


Bg^gsint  - i (Bq22+  Bq^^: 


(52c) 


It  is  clear  that  by  choosing  appropriate  combinations  of  upper  and  lower 
bounds  on  the  elements  of  B  in  Eqs .  (52),  we  can  derive  upper  and  lower 
bounds  on  \|/,  T]   ,    and  t]  . 


6.   MJMERICAL  EXAMPLES 
In  order  to  test  the  practicality  of  the  method,  we  have  applied  it 
to  the  problem  of  one  dimensional  scattering  by  an  attractive  square  well 
which  is  asymmetrically  placed  with  respect  to  the  origin,  that  is, 
V(x)  =  0  for  X  <  0  and  for  x  >  a,  and  V(x)  =  (-«  /2m)U  for  0  S  x  g  a, 
where  U  is  a  positive  constant.   (it  should  be  emphasized  that  our  method 
is  not  restricted  to  potentials  which  vanish  identically.)   Eqs.  {h)    and 
(6)  then  give 
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(55a) 


,         a  <  X  <  00 


(55b) 


where  the  range  of  x  is  now  0  to  ^.      The  trial  function  within  the  region 
of  the  square  well  was  chosen  to  be 


2P 
p=o,2,  .  .. 


u  (x) 


b  (x/a)P 
P 


2P+1 

S     b,(^/a)^ 
P=l,5,  .-.    ^ 


,  0  S  X  g  a  .      (5^) 


The  illogicalness  of  the  unnecessarily  restricted  form  of  power  series 
was  not  recognized  until  after  the  numerical  calculations  had  been  com- 
pleted, and  it  was  not  then  considered  worthwhile  to  repeat  the  calcula- 
tions with  a  less  restricted  form. 

For  all  X  >  a,  u  was  chosen  to  have  the  form  given  by  Eq.  (1^). 
u  must  of  course  be  continuous  in  slope  and  value  at  x  =  a. 

The  comparison  potentials  which  generate  non-negative  and  non-posi- 
tive matrices  A  V  were  taken  to  be 


V  (x)  =  -(^  /2m)U  1 

MA  C  Uw 


0  ^  X  ^  a 


=  0 


a  <  X  <  00 


(55) 
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and 


V  (x)  =  0  , 


0  g  X  S 


(56) 


respectively.   p  was  chosen  as 


p=p  1,     OSxSa   J   p=0 


a  <  X  < 


where  p  is  a  constant.   The  eigenphase  shifts  6  (|a)  and  5  (i-i)  corres- 
ponding to  the  lower  bound  are  just  the  phase  shifts  associated  with 
the  even  and  odd  solutions^  respectively,  of  the  scalar  differential 
equations 

^-^  +  k^s;$(x)  +  (2m/«^)nP^/(x)  =  0,    0  S  x  g  a        (57a) 
dx 


dx 


a  <  x  < 


(57b) 


while  the  eigenphase  shifts   corresponding  to  the  upper  bound  are  obtained 
from  the   same  equation  under  the   substitution 


\iP     — >     lap     +   (*  /2m)U. 
"^  o  o 


(58) 


The  phase  shifts  corresponding  to  Eq.  (57)  are  given  by 


5'    (n)   =  rt/2  -  ka  +   cot""^ 


(ka//<a)cot/<a 


(59a) 


5  '    (u)   =   -ka  +  tan 
o^ 


(ka/K:a)tanK:a 


(59b) 


2^h 


K  =(\.     +  2m|jp  /-R  J 


(60) 
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2   2 
If  we  plot  6'  ((i)  and  6'  (fj.)  versus  2in)aP  a  /■*!  ,  we  obtain  curves 

each  of  which  lies  entirely  below  the  corresponding  curve  for  the  exa.ct 

eigenphase  shift.   If  instead  of  Eq.  (6o)  we  define  K   by 


^  (k^  +  2m|ap^/«^  +  U^  ^^ 


(61) 


we  obtain  curves  which  lie  entirely  above  the  corresponding  exact  curves. 
These  pairs  of  bounding  curves  are  then  used  to  obtain  lower  bounds  on 
a  and  on  (3   in  accordance  with  the  procedure  described  in  Section  h- . 
The  particular  value  chosen  for  p  plays  no  role  whatever.   ()a  and  p 

O  «At 

appear  only  in  the  combination  |_ip.   Choosing  p  larger  would  then  give 

smaller  values  of  a  and  of  ^    ,   but  of  course  the  bounds  on  e  Be 

o       y  u*y  4*  y4*y 

would  not  be  affected. ) 

We  will  now  present  some  numerical  results  obtained  with  the  use  of 
the  Bendix  G-15  digital  computer  at  the  General  Telephone  and  Electronics 
Laboratories,  Bayside,  Long  Island.   Bounding  curves  on  S  ((a)  and  6  (|a) 

calculated  from  Eqs .  (59)  and  (6l)  are  plotted  in  Figures  3  and  k-   for 

p  op 

ka  =  0.1,  Ua  =1.0  and  in  Figures  5  and  6  for  ka  =  ^A,  Ua  =  n  /if. 

The  lower  bounding  curves  in  Figures  5  and  h   apply  also  to  ka  =  0.1, 

2 
Ua  =  0.01,  but  in  this  case  the  upper  bounding  curves  are  so  close  to 

the  lower  boxmding  curves  that  the  two  curves  cannot  be  distinguished 

from  one  another  on  the  scale  of  the  drawings.   Values  of  0  and  lower 

bounds  on  a  and  (3   consistent  with  those  curves  are  listed  in  Table  1. 
y     y 

In  each  case,  9  was  somewhat  arbitrarily  chosen  to  provide  approximately 

2 
equal  lower  boionds  on  a  and  p  .   For  each  choice  of  ka  and  of  Ua  ,  three 

y     y 

choices  were  made  for  e  ,  and  upper  and  lower  bounds  were  calculated  on 

y 

the  corresponding  values  of  e  B^e  .   The  three  choices  of  e  were 

4n  y  M.  y  im9  «  y 
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Fl gure  3 

2 

Upper  and  lower  bounds  on  5^(  [i)  for  ka  =  0.1,  Ua  =1.0.   The 

2 
lower  bound,  S  .(|^),  also  applies  to  ka  =  0.1,  Ua  =  0.01. 

The  corresponding  upper  bound,  5^^(n),  cannot  be  distinguished 

from  6   (u).   The  lower  bounds  on  a.  were  actually  chosen  from 

the  numerical  data  rather  than  from  the  curves. 
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Figure  k 

2 
Upper  and  lower  tounds  on  5^(|i)  for  ka  =  0.1,  Ua  =1.0.   The 

2 
lower  boimd,  5  A\i) ,    also  applies  to  ka  =  0.1,  Ua  =0.1. 

O  V 

The  lower  bounds  on  p  were  actually  chosen  from  the  numerical 
data  rather  than  from  the  curves.   Both  bounds  were  derived 
from  the  same  curve  and  could  have  been  chosen  equal. 
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Figure   ^ 
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Upper  and  lover  bovmds   on  5   (ti)    for  ka  =  jt/^,   Ua     =  n   /^ 
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Figure  6 


P  P 

Upper  and  lower  bounds   on  6    (|a)    for  ka  =  n/^,    Ua     =   jt    /4, 
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Table  1.   Values  of  9  and  lower  bounds  a   and  (3   for  three 
choices  of  ka  and  Ua  _,  obtained  from  Figs.  3;.  ^^  5  and.  6. 


ka 

Ua2 

9 

^91 

^Ql 

0.1 

0.01 

3-08559 

5.0 

5.0 

0.1 

1.0 

3-08359 

h.o 

h.o 

n/h 

n'/k 

2.775 

3-3 

3-5 
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The  corresponding  values  of  e  B.e   are  B^ ^ ,  B__,  and  B^^+  B__  +  2B^„. 

*  9  M  9m9       1±    dd.  11    dd.  Id 

While  it  has  been  noted  that  three  different  values  of  e   are  reqirired 

*M  9 

in  order  to  obtain  bounds  on  all  three  independent  elements  of  B  ^  no 

k*  9 

criterion  has  been  given  for  selecting  these  values.   It  is  probably 
desirable  to  bound  B  ,  and  B   separately  but  the  remaining  value  of  e 

11        did  iM  9 

could^  in  principle,  be  selected  to  provide  the  closest  bounds  on  B, „• 

However,  in  the  calculations  presented  here  no  attempt  has  been  made 

to  optimize  the  choice  of  e  . 

*u  9 

The  calculations  were  based  on  the  trial  function  given  by  Eq.  (5^) 
with  terms  up  to  x  ,  so  that  there  were  four  linear  variational  para- 
meters.  Due  to  the  two  boundary  conditions  at  x  =  a,  this  corresponds 
to  two  free  parameters .   The  free  parameters  were  varied  to  give  the 
best  bounds.   The  bounds  obtained  on  the  various  combinations  of  the 

elements  of  B  are  listed  in  Table  2.   The  bounds  are  seen  to  be  very 

2 
much  closer  in  the  case  in  which  ka  =  0.1,  Ua  =0.01  than  in  the  other 

two  cases.   This  result  is  to  be  expected,  for  considering  k  to  be  fixed, 

this  case  involves  a  less  effective  potential  thaa  either  of  the  other 

two.   Note  that  for  the  second  and  third  sets  of  values,  the  potential 

is  by  no  means  negligible  as  can  be  seen  (Table  k)    from  the  exact  values 

of  the  phase  shifts  and  the  mixing  parameters;  for  zero  potentials,  one 

would  have  t]  =  90  ,  ti  =  0  ,  and  \|;  =  0  . 
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Table  2.   Upper  and  lower  bounds  on  combinations  of  the  elements  of  the 

2 
susceptance  matrix  for  three  choices  of  ka^  Ua  and  0,  compared  with  the 

exact  values.   Two  parameters  were  used  in  the  calculations,  and  lower 

boimds   on  a     and   (3     were   determined  from  Figs.    5,    k,    5   and  6. 


quantity 


upper 
bound 


lower 
bound 


exact 
value 


ka  =   0.1,    Ua     =0.01 


B, 


11 


B, -,    +  Bp2  +   2B-,  p 


0.10803l'^■ 
-17.1726221 
-16.9786186 


0.1080085 
-17.1726509 
-16.9786503 


0 . 1080210 

■17.1726250 
-16.9786500 


ka  =   0.1,    Ua     =1.0 


B 


11 


B22 

^11  +  ^22  +  ^\2 


i+.  87112 

-12.95757 
-0.6^1^52 


^.^3705 

-15.08907 

-I.l4i68 


ij-.  71216 

-12.98585 

-0.85519 


2        2 
ka  =  ]t/4,    Ua     =   rt   A 


\l 
B 


22 

\l  +  ^22  "^   ^^12 


1.76if88 

-1.22869 

2.52112 


1.25977 
-1.51261 

1.75^95 


1.58228 

■1.25519 
2.09091 
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2    2 
Table  3-   Bounds  on  B   for  ka  =  jtA,  Ua  =  -a   /h   for  different  numbers 

of  parameters.   The  values  of  0,  a   and  p   are  as  shown  in  Table  1. 

y  1-'      y  V 


number  of 

free 

upper 

lower 

exact 

parameters 

bound 

bound 

value 

2 

1.76^^88 

1.23977 

1.58228 

4 

1.69009 

1.45289 

1.58228 

6 

1.65803 

1.49891 

1.58228 
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Table  k .      Upper  and  lower  bounds  on  the  eigenphase  shifts  and  mixing 

2 
parameters  for  three  choices  of  ka,  Ua  and  0^  compared  with  the  exact 

values.   These  quantities  were  calculated  from  the  bounds  on  the  com- 
binations of  the  elements  of  B  listed  in  Table  2,    using  Eqs .  (62). 


quantity 


upper 
bound 


lower 
bound 


exact 
value 


ka  =  0.1,  Ua 


0.01 


^ 
\ 


-0°  8'52.90"   -0°  8'53.29"   -0°  8'33.08" 


92^50 '56" 


92  50 '50 


iRr^" 


92  50 '55" 


0°  0'5^.28"    0°  0'5^.24"    0°  0'54.26' 


ka  =  0.1;,  Ua 


1.0 


-io°5^' 

-12°16' 

-11°25' 

16T°12' 

165°  r 

l66°l8' 

1°  o> 

o%r 

0°51' 

2     2 

ka  =  rtA,  Ua  =  n   /h- 


-lo°i+8' 


157°27' 


-22°  k' 


116°  y 


-15°56' 
150°24 ' 


27°53' 


1°30' 


11°3^' 
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The  effect  of  increasing  the  niomber  of  parameters  is  shown  in 

2    2 
Table  5^  which  shows  bounds  on  B   for  the  case  ka  =  n/'+j  Ua  =  :rt  /^, 

calculated  with  2,    h ,    and  6  free  parameters.   Wliile  the  bounds  improve 

as  the  number  of  free  parameters  is  increased,  the  improvement  is  not 

as  rapid  as  one  might  expect.   This  result  is  attributed  to  the  poor 

choice  of  the  form  of  the  trial  function.   The  trial  function  used  in 

calculating  upper  bounds  on  B   with  2  parameters  is  compared  with  the 

exact  function  in  Figure  7- 

Finally,  bounds  on  the  eigenphase  shifts  and  mixing  parameters 
are  listed  in  Table  k .      These  quantities  were  calculated  from  the  boionds 
on  combinations  of  the  elements  of  B  listed  in  Table  2,    using  Eqs .  (52). 
It  would  be  a  trivial  matter  to  also  obtain  bounds  on  the  elements  of  S, 
but  the  title  of  the  paper  notwithstanding,  we  have  not  bothered  to 
do  so. 

From  the  foregoing  results,  it  is  clear  that  useful  bounds  can  be 
calculated  on  the  independent  quantities  which  determine  the  scattering 
matrix  in  a  multi- channel  problem,  even  with  relatively  extreme  values 
for  the  depth  and  range  of  the  potential.   Furthermore,  since  diagonal 
comparison  potentials  can  be  used,  the  calculation  requires  no  more  than 
the  exact  solution  of  related  single- channel  problems.   Thus  it  is  anti- 
cipated that  the  method  presented  here  can  provide  useful  results  in  many 
multi-channel  problems  for  which  there  would  otherwise  be  no  practical 
approach. 
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Figure  7 


Comparison  of  trial  function  with  exact  function.   The  trial 
fimction,  which  contains  two  free  parameters,  was  used  to 
calculate  an  upper  bound  on  B^  for  the  case  in  which  ka  =  jtA, 
Ua^  =  /jh. 
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